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Abstract 

In this paper, we study adaptive finite element approximations in a 
perturbation framework, which makes use of the existing adaptive finite 
element analysis of a linear symmetric elliptic problem. We prove the 
convergence and complexity of adaptive finite element methods for a class 
of elliptic partial differential equations. For illustration, we apply the 
general approach to obtain the convergence and complexity of adaptive 
finite element methods for a nonsymmetric problem, a nonlinear problem 
as well as an unbounded coefficient eigenvalue problem. 
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1 Introduction 

The purpose of this paper is to study the convergence and complexity of adap- 
tive finite element computations for a class of elliptic partial differential equa- 
tions of second order and to apply our general approach to three problems: a 
nonsymmetric problem, a nonlinear problem, and an eigenvalue problem with 
an unbounded coefficient. One technical tool for motivating this work is the 
relationship between the general problem and a linear symmetric elliptic prob- 
lem, which is derived from some perturbation arguments (see Theorem 13.11 and 
Lemma 13.11 ) . 

Since Babuska and Vogelius [3] gave an analysis of an adaptive finite ele- 
ment method (AFEM) for linear symmetric elliptic problems in ID, there are 
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a number of work on the convergence and complexity of adaptive finite ele- 
ment methods in the literature. For instance, Dorfler [10] presented the first 
multidimensional convergence result, which has been improved and generalized 
in [51 H3 HI [T7J [THl HH1 HOI 123] ■ For a nonsymmetric problem, in particular, 
Mekchay and Nochetto [T7] imposed a quasi-orthogonality property instead of 
the Pythagoras equality to prove the convergence of AFEM while Morin, Siebrt, 
and Veeser [20| showed the convergence of error and estimator simultaneously 
with the strict error reduction and derived the convergence of the estimator 
by exploiting the (discrete) local lower but not the upper bound. To our best 
knowledge, however, there has been no any work on the complexity of AFEM 
for nonsymmetric elliptic problems in the literature. In this paper, we can get 
the convergence and optimal complexity of nonsymmetric problems from our 
general approach. For a nonlinear problem, Chen, Hoist and Xu [7 proved 
the convergence of an adaptive finite element algorithm for Poisson-Boltzmann 
equation while we are able to obtain the convergence and optimal complexity of 
AFEM for a class of nonlinear problems now. For a smooth coefficient eigenvalue 
problem, Dai, Xu, and Zhou [3] gave the convergence and optimal complexity of 
AFEM for symmetric elliptic eigenvalue problems with piecewise smooth coeffi- 
cients (see, also convergence analysis of a special case [THUS]). In this paper we 
will derive similar results for unbounded coefficient eigenvalue problems from 
our general conclusions, too. We mention that a similar perturbation approach 
was used in [5]. 

This paper is organized as follows. In section 2 we review some existing 
results on the convergence and complexity analysis of AFEM for the typical 
problem. In section 3 we generalize results to a general model problem by using 
a perturbation argument. In section 4 and section 5, we provide three typical 
applications for illustration, including theory and numerics. 

2 Adaptive FEM for a typical problem 

In this section, we review some existing results on the convergence and com- 
plexity analysis of AFEM for a boundary value problem in the literature. 

Let SI C M. d (d > 2) be a bounded polytopic domain. We shall use the 
standard notation for Sobolev spaces W s ' p (fl) and their associated norms and 
seminorms, see, e.g., [TJ [S]. For p = 2, we denote iP(Sl) = VF S ' 2 (S1) and 
Hq(Q) = {v e ii^ 1 (SI) : v |an= 0}, where v \qq= is understood in the sense 
of trace, || • || s ,q = || • || s ,2,n- Throughout this paper, we shall use C to denote 
a generic positive constant which may stand for different values at its different 
occurrences. We will also use A < B to mean that A < CB for some constant C 
that is independent of mesh parameters. All constants involved are independent 
of mesh sizes. 
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2.1 A boundary value problem 

Consider a homogeneous boundary value problem: 

f Lu := -V • (AVu) =/ in fi, 

\ 11 =0on 3!!, 1 J 

where A : £1 — >• M dx<i is piecewise Lipschitz over initial triangulation To, for 
a; G O matrix A(x) is symmetric and positive definite with smallest eigenvalue 
uniformly bounded away from 0, and / G L 2 (fl). 

Remark 2.1 The choice of homogeneous boundary condition is made for ease 
of presentation, since similar results are valid for other boundary conditions J§|/. 

The weak form of (|2.1j) reads as follows: Find u G Hq(SY) such that 

a(u, v) = (f,v) VveH'in), (2.2) 

where a(-,-) = (AV-, V-). It is seen that o(-, •) is bounded and coercive on 
i?o(fi), i.e., for any w,v £ i? 1 (f2) there exist constants < c a < C a < oo such 
that 

\a(w,v)\ < CalMknlMko and c alMli,n < a(v,v) W G H£(Sl). 

The energy norm || • || 0i q , which is equivalent to || ■ , is defined by ||iy|| a .n = 
\J a(w, w) . It is known that (|2.2[) is well-posed, that is, there exists a unique 
solution for any / G P _1 (r2). 

Let {Th} be a shape regular family of nested conforming meshes over f2: 
there exists a constant 7* such that 

where, for each r G Th, h T is the diameter of r, p r is the diameter of the biggest 
ball contained in t, and ft, = max{/i r : r G 7/j}. Let £h denote the set of interior 
sides (edges or faces) of Th- Let Sq(Q) C Hq(Q) be a family of nested finite 
element spaces consisting of continuous piecewise polynomials over Th of fixed 
degree n > 1, which vanish on <9f2. 

Define the Galerkin-projection P/j : Hq(CI) — > Sq(SY) by 

a{u - P h w, u) = Vt> G (fi)- (2.3) 
For any u G H^fl), there apparently hold: 

\\Phu\\ a ,n < ||w||o,n and lim ||u - P ft u|| a ,n = 0. 
Now we introduce the following quantity: 

Pn( h ) = sup inf ||Z/ _1 / - u||a,n, 

/eL 2 (a),||/|| 0i n=i"6Sj(n) 
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then p a (h) -> as h -> (see, e.g., [IdB])- 

A standard finite element scheme for (|2.2p is: Find ut € Sq(SI) satisfying 

a(u hl «) = (/,«) VveS#(n). (2.4) 

By definition (|2.3[) . we know that = PhU. 

By a contradiction argument, we have (c.f., e.g., [30 ) 

Lemma 2.1 As operators over Hq(Q), there holds 

lim \\K(I-P h )\\ =0 

if K, is a compact operator over Hq(Q). 



2.2 Adaptive algorithm 

Given an initial triangulation To, we shall generate a sequence of nested con- 
forming triangulations Tk using the following loop: 

SOLVE ESTIMATE MARK ->■ REFINE. 

More precisely to get 7fc+i from 7fe we first solve the discrete equation to get 
Uk on Tk ■ The error is estimated using Uk and used to mark a set of elements 
that are to be refined. Elements are refined in such a way that the triangulation 
is still shape regular and conforming. We assume that the solutions of finite- 
dimensional problems can be solved to any accuracy efficiently^ Examples of 
such optimal solvers are multigrid method or multigrid-based preconditioned 
conjugate gradient method. 

Now we review the residual type a posteriori error estimators for finite ele- 
ment solutions of (|2.ip . Let T denote the class of all conforming refinements by 
bisection of To- For Th € T and any v G Sq(CI) we define the element residual 
lZ T (v) and the jump residual J e (v) by 

Kr{v) := f-Lv = f + V-{AVv) mreTh, 
J e (v) := -AV» + ■ i/ + - AVtr ■ i/ _ := [[AVv]] e • v e on e £ 4, 

where e is the common side of elements r + and t~ with unit outward normals 
v + and v~ , respectively, and v e — v~ . Let uj e be the union of elements which 
share the side e and oj t be the union of elements sharing a side with r. 
For r € Tk, we define the local error indicator fjh(v,T) by 

fjl(v,T):=h 2 T \\iZ T (v)\\l T + J2 M-Wllo.e 

and the oscillation 6sch{v,r) by 

osc h (v, t) := hl\\H T {v) -TZr(v)\\l >r + h e \\ J e (v) - J e (v) \\l e , 

e££ h: eCdr 

x By the similar perturbation argument, indeed, it will be seen that some approximations 
to the finite-dimensional problem will be sufficient. 
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where w is the I/ 2 -projection of w E L 2 (il) to polynomials of some degree on r 
or e. 

Given a subset u C fl, we define the error estimator fjh(v,u>) and the oscil- 
lation osch(v,uj) by 

fjl(v,uj):= Vh( v > T ) and oscl(v,w):= ^ osc^r). 

For r € 7h, we also need notation 

7?£(A, r) := fc?(||divA||g i00jT + ^ 2 |l A||g )0OiUT ) 

and 

os4(A, r) := ft£(||divA - divAII^^ + /j- 2 ||A - A||§ )0Oi(1 J, 

where is the best L°°-approximation in the space of discontinuous polynomials 
of some degree. 

Given a subset wc!lwe finally set 

i]h(A,uj) :— max %(A,r) and osc^(A,w):= max osc/ l (A,r). 

We now recall the well-known upper and lower bounds for the energy error 
in terms of the residual- type estimator (see, e.g., [T71 ITUl 125] ). 

Theorem 2.1 (Global a posterior upper and lower bounds). Let u e -ffo(^) 6e 
ifte solution of \2. 2]) and Uh € Sq(SI) be the solution of \2.$ - Then there exist 
constants C\ , C2 and C3 > depending only on the shape regularity 7* , C' a and 
c a such that 

\\u-u h \\l n <C 1 fj 2 h (u h ,T h ) (2.5) 

and 

C 2 fjl(u h ,Th) < \\u - u h \\l t a + C 3 osc h (u h ,Th)- (2.6) 

We replace the subscript h by an iteration counter called k and call the 
adaptive algorithm without oscillation marking as Algorithm Do, which is 
defined as follows: 

Choose a parameter < 8 < 1 : 

1. Pick any initial mesh To, and let k = 0. 

2. Solve the system on To for the discrete solution uo- 

3. Compute the local indicators fjk- 

4. Construct ftAk C Tk by Marking Strategy Eo and parameter 9. 

5. Refine Tk to get a new conforming mesh Tk+i by Procedure REFINE. 
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6. Solve the system on Tk+i for the discrete solution itfc+i- 

7. Let k = k + 1 and go to Step 3. 

The marking strategy, which we call Marking Strategy Eq, is crucial for 
our adaptive methods. Now it can be stated by: 
Given a parameter < 9 < 1 : 

1. Construct a minimal subset Aik of Tk by selecting some elements in Tk 
such that 

f) k {u k ,Mk) > 0fjk(uk,Tk)- 

2. Mark all the elements in Aik- 

Due to [6] , the procedure REFINE here is not required to satisfy the Interior 
Node Property of [13 [19]. 

Given a fixed number b > 1, for any Tk £ T and a subset Mk C Tk of marked 
elements, 

%+i =REFINE(T fe ,M fc ) 

outputs a conforming triangulation Tk+i € T, where at least all elements of Aik 
are bisected b times. We define Rj- k ^j- k+1 :— Tk\(TkC\Tk+i) as the set of refined 
elements, thus Aik C Rr k ^T k +i- 

Lemma 2.2 (Complexity of Refine). Assume that To verifies condition (b) of 
section 4 in \2$ . For k > let {7fe}fc>o & e an D sequence of refinements of To 
where Tk+i is generated from Tk by Tk+i = REFINE(7fc, M. k ) with a subset 
MkC%. Then 

k-l 

#T fe -#To<^#VW j Vfc>l (2.7) 

3=0 

is valid, where the hidden constant depends on To and b. 

The convergence of Algorithm D is shown in [5]. 

Theorem 2.2 Let {uk}keN be a sequence of finite element solutions corre- 
sponding to a sequence of nested finite element spaces {Sp (^)}fceN produced by 
Algorithm Do- Then there exist constants 7 > and £ G (0, 1) depending only 
on the shape regularity of meshes, the data and the marking parameter 0, such 
that for any two consecutive iterates we have 

II"- «*+i|lo,n + 7Vk+x( u k+i,T k+ i) 
< i 2 (\W-Uk\\ln+lvl(uk,Tk)). 
Indeed, constant 7 has the following form 

7:= (l + 5-i)A 1% 2 (A,T )' (2 - 8) 

where t]q(A,To) ■— rff {K,To), Ai := (d+ l)Cg/c Q with Co some positive con- 
stant and constant S € (0, 1). 
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Following [51 [3] , we have a link between nonlinear approximation theory and 
the AFEM through the marking strategy as follows. 

Lemma 2.3 (Optimal Marking). Letuk & Sq(Q) anduk+i G Sq +1 (Q) be finite 
element solutions of A2.ty) over a conforming mesh Tk and its refinement Tk+i 
with marked element Aik- Suppose that they satisfy the decrease property 

\W - Uk+i\\l,n + f*6scfe +1 (u fe+ i,7fe+i) 

< /3* - u k \\l t n + i*6scl(u k ,Tk)) 

with constants % > and 0„ € (0, yi)- Then the set 1Z := RT k ^T k +i satisfies 
the following inequality 

fjk{uk,TZ) > 0fik{uk,Tk) 
mth ^ = c a (dili+2cX)i,) ' where C = A i° sc §(A,T ) and C = max(l, 

3 A general framework 

Let u E Hq(Q) satisfy 

a(u,v) + (Vu,v) = (lu,v) We#o(fi), (3.9) 

where £ : H^(fl) — > L 2 (f2) is a bounded operator and V : H^(fl) — » L 2 (fi) is a 
linear bounded operator. 

Let K : L 2 (Q) —> Hq(Q) be the operator defined by 

a(Kw 7 v) = (w,v) Vu>,v G L 2 (f2). 

Then K is a compact operator and (|3.9I) becomes as 

u + KVu = K£u. 

Let Uh € Sq(VL) be a solution of disctetization 

a(u h ,v) + (Vu h ,v) = (£ h u h ,v) Vw e Sf(n), (3.10) 

where : Sq(H) — > L 2 (f2) is some bounded operator. Note that we may view 
as a perturbation to for which we assume that there exists Ki(h) <G (0, 1) 
such that 

\\K(£u-£ h u h )\\ a ,n = 0(K 1 (h))\\u-u h \\ a ,n, (3.11) 

where Ki(h) — > as h — > 0. 

Note that (|3.10p can be written as 

u h + P h KVu h = P h K£ h u h , 

where Ph is defined by (|2.3|) . We have for w h = K£hUh — KVuh that 

u h = P h w h . (3.12) 
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Theorem 3.1 There exists n(h) £ (0, 1) such that n(h) — > as h — > and 

\\u - u h \\ a . n = \\w h - P h w h \\ a>n + 0{n{h))\\u - u h \\ a , n . (3.13) 

Proof. By definition, we have 

u-w h = K£u - KVu - (K£ h u h - KVu h ) = K(£u - £ h u h ) + KV{u h - u). 

Let K 2 (h) = \\KV(I-P h )\\. Since KV : H^(tt) H£(fL) is compact, we get 
from Lemma H3] that K 2 (h) — > as h — » 0. Note that 

KV(u h - u) = KV(I - P h )(u h - u), 

we obtain 

\\KV(u h - u )|| tti0 = 0{k 2 {K))\\u - Ufc || Bi n. (3.14) 
Set k(K) = K\(h) + K 2 (h), we have that k(K) — > as h — > and 

\\u - w h \\ at u < Ck(K)\\u - u h \\ a ,n. (3.15) 

Since (|3.12p implies 

u - u h = w h - PhW h + u — w'\ 

we get (|3.13l) from p. 151) . This completes the proof. 

Theorem 13.11 sets up a relationship between the error estimates of finite el- 
ement approximations of the general problem and the associated typical finite 
element boundary value solutions, from which various a posteriori error estima- 
tors for the general problem can be easily obtained since the a posteriori error 
estimators for the typical boundary value problem have been well-constructed. 
In fact, Theorem 13.11 implies that up to the high order term, the error of the 
general problem is equivalent to that of the typical problem with £h,Uh — Vuh 
as a source term. However, the high order term can not be estimated easily in 
the analysis of convergence and optimal complexity of AFEM for the general 
problem, for instance, for a nonsymmetric problem, a nonlinear problem and an 
unbounded coefficient eigenvalue problem. 

3.1 Adaptive algorithm 

Following the element residual lZ T (u,h) and the jump residual J e (uh) for (|2.4p . 
we define the element residual lZ T (uh) and the jump residual J e {uh) for (|3.10l) 
as follows: 

Tlriuh) := £ h u h - Vu h - Lu h = £ h u h - Vu h + V • (AVti/,) in r e %, 
J e (u h ) := -AV4 • v + - AVu ft • v~ := [[AV%]] e ■ Ve on e e £ h - 

For r £ 7~h, we define the local error indicator rjh(uh, t) by 

r)l(u h ,T):=h 2 T \\K T (u h )\\l T + £ MI-MOIIo.e 
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and the oscillation osch{uh,T) by 

os(? h {u h ,T) := h 2 T \\U T {u h ) - H T (uh)\\l yT + ^ h A\ J e{uh) - Je.{uh)\\Q, e , 

where e , v + and v~ are defined as those in section [51 

Given a subset u c SI, we define the error estimator rjh(uh,w) by 

rfc(u h ,u):= 4(uh,r) (3.16) 

and the oscillation osc^(m^,w) by 

osc 2 h {u h ,u) := ^ os4(u ft ,r). (3.17) 

t£Ti,,tCu 

Let /io S (0,1) be the mesh size of the initial mesh 7o and define 

R,(ho) := sup n{h). 

he{0,h o ] 

Obviously, k(ho) 1 if ho <C 1. 

To analyze the convergence and complexity of finite element approximations, 
we need to establish some relationship between the two level approximations. 
We use Th to denote a coarse mesh and Th to denote a refined mesh of Th- 
Recall that w h — K(£hUh — Vuh) and w H = K(Ihuh — Vuh)- 

Lemma 3.1 Let h,H £ (0, ho], then 

\\u- u h \\ a ,n = \\w H ~ P h w H \\a,n + O(R(h )) Q\u - u h \\ a>n + \\u- u H \\ a ,n) ,(3.18) 

r)h(uh,%) = ijh(PhW H ,%) + O{h{h )) {\\u- u h \\ a> n + \\u- u H \\ a ,n) , (3.19) 
and 

osc h (u h ,%) = osc h (P h w H ,%) + O(k(h )) (\\u - u h \\ a ,n + \\u - u H \\ a ,n) (3.20) 

Proof. First, we prove (|3.18p . It follows that 

\\P h (w h -w H )+u-w H \\ a , n < \\w h -w H \\ a ,ci + \\u-w H \\ a ,n 

< \\u - w H \\ a:n + \\u - w h \\ a ,n, 

which together with (|3.15|) implies 

\\P h {w h - w H ) + w H - u\\ a . n < k(H)\\u - u H \\ a ,n + n(h)\\u - u h \\ a<n . 
Namely, 

\\P h (w h - w H ) + w H - u||„,o < R(h )(\\u - u H \\a,a + \\u- u h \\ a ,a). (3.21) 



9 



Observing that identity (13. 12)) leads to 

u-u h = w H - P h w H + P h {w H -w h ) + u- w H , 

we then obtain (|37T8|) from (pHTjl . 

Next, we turn to prove (|3.20[) . Due to Lw h = ItUh — Vut and Lw = 
^huh — Vuh , we know that w h — w H is the solution of typical boundary value 
problem with thUh — £rUh + Vuh — Vuh as a source term. Since 

H T (P h (w h - w H )) = t h u h - l H u H + Vu H - Vu h - L(P h (w h - w H )), 

we have 

5Ic 2 h (P h {w h -w H ),T h ) = ^ c I(E,t) 

reT h 



^2 (h 2 T \\1Z T (E)-n T (E)\\l T + J2 K\\J e (E) - J e (E)\\l e 

reTh e££ h ,e(ldT 



< J2 h r\\^r(E) +LE- (K T (E) + LE)\\l 



reT h 



+ J2( h r\\LE-LE\\l r + h e \\J E (E)-J e (E)\\l e ), (3.22) 

where E = Ph(w h — w H ). Following the proof of Proposition 3.3 in [BJ, we see 

that 

J2 (h 2 T \\LE -LE\\l T + J2 h e \\J e (E) - J e (E)\\l e ) 

can be bounded by 

C> W2(A,r)||P,(^-^)||^< OSC 2 JA,r,)||P,(^- W ff )||^. 

reT h 

Hence using the fact osCh(A.Th) < osc (A,Tq), we obtain 



J2 (h 2 T \\LE-LE\\l T + J2 h e \\ME)- Je(E)\\le) 

TeTh ee£ h ,eC9r 

< oscP (A,T )\\P h (w h -w H )\\l n . (3.23) 
Using the inverse inequality, the bounded property of V and p. lip , we get 

( %\\KtW + LE - (n T (E) + LE)\\l T ) 1/2 

Ten 

< \\hAe h u h - £ h uh + Vu H -Vu h )\\l T ) 1/2 
Ten 

< \\K(£ h u h - £HU H )\\a,n + h\\u H - u h \\ a ,n 

< \\K(£ h u h - lu)\\ a . n + \\K(i H u H - lu)\\ at a 

+h\\u- U H \\a,Q + h\\u - Uh\\a,Q 

< k(h ) (\\u - u h \\ a ,n + \\u - u H \\ a ,n) ■ (3.24) 
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Note that 

\\P h (w h -w H )\\ a> n<\\w h -w H \\ at n 

which together with (|3.15l) implies 

\\P h (w h -w H )\\ a , n < K(M(ll«-«fclla,n + ||u-«H||a^). (3.25) 

Combing $£Z5§ , (|3TM1) and ([535]) , we conclude that 

osc h {P h (w h - w H ),Th) < H(h ) (]\u- u h \\ a ,n + \\u - u H \\ a ,n) ■ (3.26) 

Due to Uh = PhW H + Ph(w h — w H ), we obtain from the definition of oscillation 
that 

oTc h {P h w h ,T h ) < oTc h {P h w H ',%) + dTc h (P h (w h - w H ),T h ). (3.27) 

Hence from osch(uh,Th) = osch(uh,Th), (I3.26|) and (|3.2T[) . we arrive at (|3.20l) . 
Finally, we prove (p7l9)l . By (HTSJ) and (|3~26|) . we have 

Vh(P h (w h -w H ),T h ) 

< \\(w h -w 11 )- P h (w h - w H )\\ a , n + ™c h {P h {w h -w H ),T h ) 

< \\u - w h \\ a: n + \\u - w H \\ a ^ + k(h ) (\\u - u h \\ a ,n + \\u - u H \\ a ,n) 

< R(ho) (\\u - u h \\ a ,Q, + \\u - u H \\ a ,n) • (3.28) 

From (j3T^51) and the fact that 

Vh{P h w h 7 T h ) = v h {P h w H + P h {w h - w H ),T h ), 

we obtain 

fj h (P h w h ,Th) = i) h {P h w H ,%) + O(k(h )) (\\u - u h \\ a ,n + \\u - u H \\ a ,n) , 

which is nothing but (|3.19[) since f}h{Ph,w h ,Th) = Vh{uh,Th)- 

Theorem 3.2 Let ho -C 1 and h € (0, ho]. There exist constants C\^Ci and 
C3, which only depend on the shape regularity constant 7*, C a and c a such that 

\\u-u h f a<n <Cxri 2 h (u h ,T h ) (3.29) 

and 

C 2 Vh(uh,T h ) < \\u-u h \\l n + C 3 osc 2 h (u h ,T h ). (3.30) 
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Proof. Recall that Lw h — IhUh — Vut- From (I2.5[) and (|2.6|) we have 

- PhW h \\l,n < Cxr,l{P h w h ,T h ) (3.31) 

and 

C 2 f h (P h w h ,T h ) < \\w h -P h w h \\l n + C 3 oTc 2 h (P h w h ,T h ). (3.32) 

Thus we obtain and (pHO]) from f3~T2j) . (f3A3|) . (j3~3T]) and ([3~32"1) . In 

particular, we may choose Ci, C2 and C3 satisfying 

Ci -^(l + C^o)) 2 , C 2 = C 2 (1-Ck(h )) 2 , C 3 = C 3 (l-Ck(h )) 2 (3.33) 



Remark 3.1 The requirement ho <C 1 is somehow reasonable for finite element 
approximations of &3.9\) . We can refer to [111 for the initial mesh size require- 
ment in adaptive finite element computations for nonsymmetirc boundary value 
problems. 

Now we address step MARK of solving (|3.10p in detail, which we call 
Marking Strategy E. Similar to Marking Strategy Eq for (|2.4p . we de- 
fine Marking Strategy E for (|3 . 10[) to enforce error reduction as follows: 

Given a parameter < 9 < 1: 

1. Construct a minimal subset Aik of Tk by selecting some elements in Tk 
such that 

r/ k (u k ,Mk) > On k (u k ,Tk)- 

2. Mark all the elements in Ai k - 

The adaptive algorithm of solving (|3.10p . which we call Algorithm D, is 
nothing but Algorithm Do when Marking Strategy Eo is replaced by Mark- 
ing Strategy E. 

3.2 Convergence 

We now prove that Algorithm D of (|3.10j) is a contraction with respect to the 
sum of the energy error plus the scaled error estimator. 

Theorem 3.3 Let 8 G (0, 1) and {uk}kef) oe a sequence of finite element so- 
lutions corresponding to a sequence of nested finite element spaces {Sq (^l)}k^ 
produced by Algorithm D. Then there exist constants 7 > and £ € (0,1) 
depending only on the shape regularity constant 7*, C a , c a and the marking 
parameter 9 such that 

\\ u ~ u k+i\\lv +7»?fc+i(«fc+i,7fc+i) 
< e 2 (h-"fc|la,o+7^K,T fc )). (3.34) 
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Here, 

with C4 a positive constant, provided /io<l. 

Proof. For convenience, we use uh, uh to denote Uk+i and respectively. 
Thus we only need to prove that for Uh and uh, there holds, 

\\ u ~ Uh\\l,n +7Vh( u h,Th) < £, 2 {\\u-UH\\l,n+VlH( u H,TH))- 

We conclude from Theorem 12.21 w h — K(lhUh — Vuh) and w H — K(£huh — 
Vuh) that there exist constants 7 > and £ <G (0, 1) satisfying 

\\w H - p h w H \\i n + ^i{p h w H ,r h ) 

< i 2 (\\w H P H w H \\l n +jfj 2 H (P H w H ,Th)). 

Hence use the fact that uh — Phw H , we obtain 

\\w H ^ P h w H \\l n + ^ 2 h (P h w H ,T h ) 
< i 2 {\\w H - u H \\ln +1Vh(vh,Th)). (3.36) 

By p,18p and (|3.19p . there exists a constant C > such that 

II" - u h\\l,n+Tlh(uh,T h ) 
< (1 + §i)\\w H - P h w H \\l n + (1 + Si^fjliPhw" ,Th) 
+C(1 + S^)k 2 (h Q )(\\u - UhH^n + ||u - u H \\l n ) 
+6(1 + S^)~K 2 (h a M\u- u h \\l t{l + \\u- u H \\l t n), 

where the Young's inequality is used and Si € (0, 1) satisfies 

(1 + <5i)£ 2 <1. (3.37) 

It thus follows from (|3.36p . (|3.15p , and identity t]h{Phw H ,Th) = Vh(uh,Th) 
that there exists a positive constant C* depending on C and 7 such that 

\\ u - u h\\l,a + lVh( u h,Th) 

< (1 + S 1 )i 1 {\\w H - u H f a<n + jvh(u h , T h )) 
+C*S^ 1 k 2 (h ){\\u - u h \\ 2 a U + \\u - u H \\l t a) 

< (1 + S^ 2 ((1 + Ck(h Q )) 2 \\u - u H \\l n + ^ 2 H (u H ,T H )) 
+C*S^k 2 {h ) (\\u - uh\\ 2 at a + \\u- u H \\l,a) ■ 
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Hence, if ho <C 1, then there exists a positive constant C 4 depending on C* and 
C such that 

\W - u h \\l >Q + tyriliu^Th) 
< (1 + <5i)f (||u - u H ||2,n + JV 2 h(uh,Th)) 

+C 4 k{h )\\u - u H \\ 2 a n + C 4 Si 1 K 2 (h )\\u - u h \\ 2 a fl . 

Consequently, 

(1 - C 4 5i l K 2 (ho)) \\u - u h \\ 2 an + JVh( u h,%) 
< ((1 + Sx)i 2 + C 4 k(h )) \\u - ua\\lfl + (1 + h)?ir? H {u H ,T H ), 

that is 

\\ u ~ u h\\l,n + „ Zi~ 2lu z vU u h,T h ) 

1 — C 4 o 1 K^{no) 



(l + S 1 )e + C i R(h ) 2 (l+fr)£ 2 7 

l-C^K^/lo) l-C^R 2 ^)' 



2 



Since /iq -C 1 implies f(/io) <C 1, we have that the constant £ defined by 

1/2 

/ ( I -!- A i< - !- < ',;,•!/),, 1 \ 



(i + s 4 )e + c 4 k(h ) 



1 - C 4 «f 1 K 2 (ft ) / 

satisfying £ g (0, 1) if /io -C 1. Therefore, 

ll u - u h\\l + - — „ Zi~ 9ru , vl(uh, %) 

< e (\\u-u H \\i n + -y H (u H ,r H ) 

\ (1 + S 4 )e + C 4 k(Iiq) 

Finally, we arrive at (|3.34l) by using the fact that 

(1 + 



(1 + s^e + c 4 k(h Q ) 

This completes the proof. 



< 7- 



3.3 Complexity 

We shall study the complexity in a class of functions defined by 

A* := {v £ H^Q) : |w| s>7 < 00}, 
where 7 > is some constant, 

M s , 7 = sup£ inf (#Tk-#T o y 
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and 71- C To means 7fc is a refinement of To- It is seen from the definition 
that, for all 7 > 0, At — A\. For simplicity, here and hereafter, we use A s 
to stand for A\, and use \v\ s to denote |v| s ,7- So A s is the class of functions 
that can be approximated within a given tolerance e by continuous piecewise 
polynomial functions over a partition Tk with number of degrees of freedom 
#% - #T < e-V»\v\V s . 

In order to give the proof of the complexity of Algorithm D for solving 
(|3.10p . we need some preparations. Recall that associated with Uk, the solution 
of (|3.10p in each mesh Tk, w k — K(£kUk — Vuk) satisfies 

a(w k ,v) = (£ k u k -Vu k ,v) Vt>efl_(fi). (3.38) 

Using the similar procedure as in the proof of Theorem 13.31 we have 



Lemma 3.2 Let Uk and Uk+i be discrete solutions of (3.10\) over a conforming 
mesh Tk and its refinement Tk+i with marked set Ai k . Suppose that they satisfy 
the following property 



\W - Mfc+i IL.fi + l*osc k+1 {u k+l ,Tk+i) 
< 0ti\\u - Uk\\l,Q + j*osc 2 k (u k ,T k )) , 

where 7* and /3* are some positive constants. Then for problem US. 38\) , we have 

\\w k -P k+1 w k \\l n +^oYc 2 k+1 (P k+1 w k ,T k+ i) 

f3*\\\w k ~P k w k \\l n + i*olfc 2 k (P k w k ,T k )) 



< 



with 



l-C 5 S^k 2 (h ) J ' l-C 5 S^k 2 (h y 

where C5 is some positive constant and Si e (0, 1) is some constant as in the 
proof of Theorem [ 



P* ■= — ; ^ r _i-->M \ > 7* : = ; „ .-i-->M x ' I 3 - 39 ) 



Corollary 3.1 Let u k and u k +\ be as those in Lemma WM . Suppose that they 
satisfy the decrease property 

II"- u k+i\\l,n + l*oscl +1 (uk+i,Tk+i) 
< Pl(\\u - Wfc|lo,n + l*osc 2 k (u k ,T k )) 

with constants 7* > and (3* G (0, Then the set 1Z :— RT k ^T k +i satisfies 

the following inequality 

■q k (uk,Tl) > 0t]k(uk,Tk) 
with 9 2 = - ,^ 2 ^~ 2 ^*^ \ - , and Co = max(l, &). where and 7* are defined 

Oo(Ci + (l+2GGi)7») 7* 
_____ _ 2 

in VS. 39\) with 5\ being chosen such that /3* G (0, _). 
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Proof. It is a direct consequence of combining u k — P k w k with Lemma [ 
and Lemma 13.21 

The key to relate the best mesh with AFEM triangulations is the fact that 
procedure MARK selects the marked set M. k with minimal cardinality. 

Lemma 3.3 (Cardinality of Let u G A s , 7k be a conforming partition 

obtained from To, and satisfies € (0, c 3 (Ci+{i+2CCi)-y) )• Then the following 
estimate is valid: 

#M k <(\\u-u k \\l n+1 o S cl(u k ,T k )y 1/2s \u\l/ s , (3.40) 



where the hidden constant depends on the discrepancy between 9 and 



C 2 i 



C 3 (Ci + (l+2CCib) 



with C defined in Lemma \2. 
Proof. Let a, a\ £ (0, 1) satisfy ct\ £ (0, a) and 

a < ^27 /, _ 2n 

< C 3 (d + (l + 2CCi) 7 ) 1 " ; ' 

Choose Si € (0,1) to satisfy (|3.37[) and 

{l + S^al < a 2 , (3.41) 



which implies 



Set 



(l + «i)a?<l. (3.42) 

«i(||w - «fe||a,n +Josc 2 k (u k ,T k )) 1/2 



V2 

and let 7^ be a refinement of 7o with minimal degrees of freedom satisfying 

\\u-u s \\l n + ( 1+ l)osc 2 (u e ,Te)<£ 2 . (3.43) 

It follows from the definition of A s that 

#T e - #T < e- x/ -|«|y-. 

Let 7i = 7^ © 7fc be the smallest common refinement of T k and 7^ . Note that 
w e = K(l £ u £ — Vu £ ) satisfies 

Lw £ = l £ u £ — Vu £ , 

we get from the definition of oscillation and Young's inequality that 

oTcl{P,w e , T ) < 2oTcl(P £ w e ,T) +2C 2 osc 2 (A,t)\\P £ w £ - P^lll^ VreT,, 

which together with the monotonicity property osc*(A,7^) < oscq (A, 7o) yields 



o^ct{P,w e ,%) < 2oTcl{P £ w £ ,%)+2C\\P £ w £ -P*w £ \\l 
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where C = Aioscq(A, 7q). Due to the orthogonality 



we arrive at 



P*™ S \\ln = IK - PeW e \\ln ~ \\P*™ £ ~ PeW £ \\la, 



IK - P*w e \\l,n + ^os&P.vf, %) 

< \\W £ - P e W e \\l n + ^0 S C 2 £ (P e W S ,Te). 



Since (|2.8[) implies 7 < tjW, we obtain that 



K - P*™ £ \\i,n + l5$K(P*w e ,%) 
1 

C~ 

< \\W £ - P e W S \\l u + ( 7 + Cr)0SC 2 e (P e W S ,Te) 



< \\w e -P e W e \\l u + -0s4(P £ W S ,Te) 



with a — ^—7 G (0, 1). Applying the similar argument in the proof of Theorem 
3. 3 1 when (|3.19[) is replaced by (|3.20p . we then get 

ll« -«* llo.fi +Joscl(u*,%) 

< «o (Ik - Mla.fi + (7 + v)osc 2 (P E w £ , %)) 

< a 2 (\\u-u e \\l n + fr + l)osc 2 e {P s w s ,T s )), (3.44) 



where 



(l + fr) + C 4 «(ho) 



i-c 4 5r^ 2 (^o) 

and C4 is the constant appearing in the proof of Theorem 13.31 Thus, by (|3.43|) 
and (|3l4|) . it follows 

\\ u - u *\\l,n + loscl(u*,%) < a 2 (\\u - u k \\ 2 an + josc 2 k {u k , %)) 

with a = -^aocti. In view of (13.42[) . we have a 2 € (0, \) when ho <C 1. Let 
1Z := Rj- k - y % t) by Corollary 13. 11 we have that % satisfies 

i] k (u k ,K) > 97] k (u kl Tk), 
Where ° 2 = cJ$wcUy 7= 1 _ c .^ sa(ho) ,go = max(l,^), an d 



2 _ (l + gi)a 2 + C 5 K(/t ) 
1 - C 5 S^k 2 (h ) 



a 



It follows from the definition of 7 (see p.35|l ) and 7 (see (|2.8p ) that 7 < 1 and 
hence (5q = S?-. Since /iq < 1, we obtain that 7 > 7 and & 6 (0, -^a) from 
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(I3.4ip . It is easy to see from p. 331) and 7 > 7 that 



C 2 {l-2a 1 ) ^ G 



fyiCx + (1 + 2C<5 1 ) 7 ) <7g(^ + i + 2CCi) 

c 2 



(l-Ck(ho)) 2 n 2 



c 3 1 Ci 1 -1 1 2c 1 

(1-C£(M) 2V 7((1+Ck(M) 2 ) ^ (l+CS(h„)) 2; 



> &2 ,.. _ 2 \ C27 /.. _ 2\ . a 

- c 3 (^ + (l + 2CC 1 )) 1 a ' C 3 (d + (l + 2Cd) 7 ) 1 " J 
when /iq <C 1. Thus 

< (^r 1/s (h - u k \\in + iosct(u k ,%)y 1/2s m;/«, 

which is the desired estimate p.40|) with an explicit dependence on the discrep- 
ancy between 9 and c 3 (Ci+(i+2CCi)y) v * a ai " This completes the proof. 
As a consequence, we obtain the optimal complexity as follows. 

Theorem 3.4 Let u € -4 s and {ufe}fcgN be a sequence of finite element solu- 
tions corresponding to a sequence of nested finite element spaces {Sq (£!)}fc 6 N 
produced by Algorithm D. Then 

\\u-u k \\l n +-yoscl(u k ,T k ) < {#% - #%)- 2s \u\l 

where the hidden constant depends on the exact solution u and the discrepancy 
between 9 and 



c 2 i 



C 3 (Ci + (l+2CCi) 7 )- 

Proof. It follows from $FT$ and (pHUl) that 



#r fc - #r < # M 3 

3=0 

fe-1 



< J2(W u - u 3\\ln+losc^u ] ,T J ))- 1/2s \u\l/ s . 

Note that (|3~3U1) implies 

\\u-v>j\\l i n+'YVj(uj,Tj) < C(\\u-u j \\l iS1 + 'yosc§(<Uj,7j)), 
where C = max(l + ^- , p 3 - ) . It then turns out 



c 2 ' c 2 < 
fe-i 



#r fc -#T < ^(\\u-uA\l# + V$(Mi,Ti)) ' Ml /a 

3=0 
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Due to (|3.34j) . we obtain for < j < k that 

\W-u k \\l >a + Ti(u k ,T k ) < e {k ~ 3] (l|w-^lla,n + 7^K-,7})) ■ 
Consequently, 

fc-i 

3=0 

< \u\\' s {\\u-u k \\l n+1 rg{u k ,T k ))- y2 \ 

the last inequality holds because of the fact £ < 1. 
Since osc k (u k ,Tk) < r) k (u k ,T k ), we arrive at 

#% - #r < (n«- Uk n^ + josct(u k ,T k )y 1/2s ny s . 

This completes the proof. 

4 Applications 

In this section, we provide three typical examples to show that our general 
theory is quite useful. 

4.1 A nonsymmetric problem 

The first example is a nonsymmetric elliptic partial differential equation of sec- 
ond order. We consider the following problem: Find u £ Hq(H) such that 

f -V • (AVu) + b • Vu + cu = f in fl, 

\ u = o on an, [ ' 

where n C M. d {d > 2) is a bounded ploytopic domain, A : ft — > M. dxd is piecewise 
Lipschitz over initial triangulation 7o, for x £ £1 matrix A(x) is symmetric 
and positive definite with smallest eigenvalue uniformly bounded away from 0, 
b € [L°°(n)] !i is divergence free , c £ L°°(n), and / e L 2 {fl) . 

A finite element discretization of (|4.45p reads: Find Uh € Sq(Q) satisfying 

{AVu h , V«) + (b • V« h ,«) + (cu h ,«) - (/, i>) V« e Sf(fi). (4.46) 

It is seen that (|4.46p is a special case of (|3.10j) . in which Vu := b • Vit + cu and 
iu = IhUh = /• Consequently, = 0, w h = K(f — Vuh) and 

u - w h = KV(u h -u) = KV(I - P h ){u h - u). 

Obviously, V : Hq(&Y) — > L 2 (f2) is a linear bounded operator and -ftTV^ is a 
compact operator over Hq(Q). We have the conclusion of Theorem 13. II 
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In this application, the element residual and jump residual become 
U T (u h ) := f - b • Vufc - cuh + V • (AViifc) in r € %, 

Je(uh) ■= [[AVuh]]e-Ve on e £ Eh 

while the corresponding error estimator r)h(V'hj Th) and the oscillation osCh(uh, Th) 
are defined by (|3.16p and (|3.17|) , respectively. Thus Theorem 13.31 and Theorem 
[33] ensure the convergence and optimal complexity of AFEM for nonsymmetric 
problem (|4.45|) . 

4.2 A nonlinear problem 

In this subsection, we derive the convergence and optimal complexity of AFEM 
for a nonlinear problem from our general theory. 

Consider the following nonlinear problem: Find u £ H^fi) such that 

f Cu := -Au + f(x,u) = in ft, UAT) 
\ u = on <9f2, ' 

where f(x,y) is a smooth function on R 3 xR 1 , 

For convenience, we shall drop the dependence of variable x in /(x, u) in the 
following exposition. We assume that u £ H^(fl) D H 1+s (f2) for some s £ (0, 1]. 
For any w £ Hq(CI) n i7 1+;5 (f2), the linearized operator £^ at iu (namely, the 
Frechet derivative of C at uu) is then given by 

C' w = -A + f'(w). 

We assume that C' w : i?o(0) — ^ H~ l {£l) is an isomorphism. As a result, u £ 
i?o(f2) ^ ff 1+s (^) must be an isolated solution of (I4.47|) . The associated finite 
element scheme for (|4.47l) reads: Find Uh £ Sq(&) satisfying 

(Vu h ,Vt;) + (/(u h ),i;) = £ S#(Q) . (4.48) 

Let o(-,-) = (V-,V-), if = (-A)" 1 : L 2 {Q) -> H^{Q), V = and 4» = -/(to) 
for any u; £ Sft(fi), then (gl5| becomes ([3TTD]) . 

As usual, to analyze the finite element approximation of nonlinear problem 
(|4.48p . we require mesh Th to satisfy that there exists <; > 1 such that (c.f. [29]) 

/i ? < h(x) xeQ, 

where h(x) is the diameter h T of the element r containing x. We consider 
the case of that Sq(SY) is the conforming piecewise linear finite clement space 
associated with Th- We assume that <, < 2s. Thus we can choose p £ (3, 6<^/(3<j — 
2s)] and obtain from Theorem 3.1 and Theorem 3.2 of that 

Lemma 4.1 If h <C 1, i/ien 

||w - «h||i,n + ^ s ||uh||o,oo : a < ft- 3 
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and 

\\u - u/Jlo.o < r(h)\\u - u h \\i,n, 

where r(h) — > as /i — > 0. 

Now we shall show that Theorem 3.1 is applicable for (|4.47[) . Since K is 
monotone and f{x,y) is smooth, we have from Lemma |4. II that 

\\K(f(u) - /K))||„,n < \\K(u - u h )\\ a>n 
< \\u - v,h\\o,n < r(ft)||it - Mh|| ,n. 

Therefore we have (|3.13p when we choose Ki(h) = r(h) and K2(h) = 0. 
In this application, the element residual and jump residual become: 

U T (u h ) := -f(u h ) + Au h in t G %, 
J e (u h ) := -Vu~l -v + -Vu^ -v~ :=[^u h ]] e -v e on e G £ h 

and the corresponding error estimator rjh{uh, Th) and the oscillation osCh{iih, Th) 
are defined by Q3.16P and p,17[) . respectively. Then Theorem 13. 31 and Theorem 
I3.4l ensure the convergence and optimal complexity of AFEM for nonlinear prob- 
lem (|PT|) . 



4.3 An unbounded coefficient problem 

Finally, we investigate a nonlinear eigenvalue problem, of which a coefficient is 
unbounded. It is known that electronic structure computations require solving 
the following Kohn-Sham equations [4j [14l [16] 

~A- V r^S + I T^Kdy + V xc (p)]u i = X i u i in R 3 , (4.49) 

1 ~t x \ x ~ r j\ Jm 3 \ x — y\ l 

where N atom is the total number of atoms in the system, Zj is the valance 
charge of this ion (nucleus plus core electrons), rj is the position of the j-th 
atom (j = l,... ,N atom ), 

N acc 
i=l 

with Ui the i-th smallest eigenfunction, Ci the number of electrons on the i-th 
orbit, and N occ the total number of the occupied orbits. The central computa- 
tion in solving the Kohn-Sham equation is the repeated solution of the following 
eigenvalue problem: Find (A, u) eKx H^fl) such that 

-^Au + Vu — Xu in 57, . 
ll«llo,n =1, 1 j 
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where CI is a bounded domain in I 3 , V = V ne + Vo is the so-called effective 
potential. Here, Vq £ L°°(fl) and 

N a tom y 

Vne(x)=- £ 

\ X - r 3\ 

A finite element discretization of (I4.50P reads: Find (Xh,Uh) £ K X Sq(CI) 
such that 

^(Vu/j, Vu) + (Vu h , v) — \h(uh, v ) Vv£S£(n). (4.51) 
Let i h : Sft(fl) -> L 2 (^) be defined by 

4,« = A h « v« £ s£(n), 

then (|4~5Tj) is a special case of (|3T0)) when a(-, ■) = |(V-, V-) and if = i(-A)" 1 : 
L 2 (Cl) Tig 1 (ft). 

Using the uncertainty principle lemma (see, e.g., [2 5) ) 

/ |V W | 2 Vw£Gr(R 3 ) 

and the fact that C£°(£l) is dense in Hq(CI), we obtain 

/^<4/|VH 2 V«Eitf(n). 

Then for any w £ i?g(f2), we have 

\\V ne w + V Q w\\ . n < C\\w\\i,n, 

namely, V is a bounded operator over Hq(CI). Thus i^U is a compact operator 
over £#(n). 

We consider the case of that (A, u) £ M x Hq (CI) is some simple eigenpair of 
(|4~50)) with ||u||o,n = !• Note that for &> := Aw Vu £ i?o(0), there holds 

- 4«h) = AAT(ti - u h ) + (A — \ h )Ku h . 

So if (Ah,Uh) £ K x Sq(CI) is the associated finite element eigenpair of (|4.51[) 
with HufrHo.n — 1 that satisfy 

|| u - Wfc||o,n + |A - A,j| < ki(/i)||u - Uft,|| a> n, 

we then have (c.f. [9]) 

||if(^u- 4uh)|| ,n = 0(Kx(h))\\u - u h \\ a ,n, 

where K\(h) := p a (h) + \\u— Uh\\ a ,n satisfying K\(h) — > as h — > 0. 
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In this application, the element residual and jump residual become: 

TZ T (u h ) := X h u h - Vu h + ^Au h in t e T h , 

J e (u h ) := [[-V«/,]] e ■ v e on e G £/j 

and the corresponding error estimator rjh{uh, 7/j) and the oscillation osch(uh,Th) 
are defined by (|3.16|) and p.!7|) . respectively. Then Theorem 13. 31 and Theorem 
13.41 ensure the convergence and optimal complexity of AFEM for unbounded 
coefficient problem (|4.50l) (c.f. 0). 



5 Numerical examples 

In this section we will report some numerical results to illustrate our theory. 
Our numerical results were carried out on LSSC-II in the State Key Laboratory 
of Scientific and Engineering Computing, Chinese Academy of Sciences, and our 
codes were based on the toolbox PHG of the State Key Laboratory of Scientific 
and Engineering Computing, Chinese Academy of Sciences. 

Example 1. We consider (|4.45j) when the homogenous Dirichlet boundary 
condition is replaced by u = g on dil and 57 ~ (0, l) 3 with the isotropic diffusion 
coefficient A = eJ, e = 10~ 2 , convection velocity b = (2,3,4), and c = (c.f. 
[15] for a 2D case and Remark l2.ip . The exact solution is given by 




For small e > the solution has the typical layer behavior in the neighbourhood 
of x = 1, y = 1, z = 1, respectively. The Dirichlet boundary condition g(x, y, z) 
on dQ is given by 

, , J x = 1 or y = 1 or z = 1, 
g\ x ->y, z ) — <y u ( x ^y >z ^ x — or y = or z = 0. 

Some adaptivcly refined meshes are displayed in Fig. 15.11 and Fig. 15.21 Our 
numerical results are presented in Fig. 15.31 and Fig. 15.41 It is shown from Fig. 
15.41 that \\u — Uh\\i is proportional to the a posteriori error estimators, which 
indicates the efficiency of the a posteriori error estimators given in section 14.11 
Besides, it is also seen from Fig. 15.31 and Fig. 15.41 that, by using linear finite 
elements and quadratic finite elements, the convergence curves of errors are 
approximately parallel to the line with slope —1/3 and the line with slope —2/3, 
respectively. These mean that the approximation error of the exact solution has 
optimal convergence rate, which coincides with our theory in section 13.21 

Example 2. Consider the following nonlinear problem: 

f -Au + u 3 = f in fl, 
1 u — on dfl, 
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Z=0.0 



Z=0.0 



Figure 5.1: The cross-section of 
an adaptive mesh of Example 1 
using linear finite elements 



Figure 5.2: The cross-section of 
an adaptive mesh of Example 1 
using quadratic finite elements 



— — the a posteriori estimator 
a lire with slope -1/3 





lumber ol degrees of freedom 



number df degrees ot freeddm 



Figure 5.3: The convergence 
curves of Example 1 using lin- 
ear finite elements 



Figure 5.4: The convergence 
curves of Example 1 using 
quadratic finite elements 
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Figure 5.5: The cross-section of 
an adaptive mesh of Example 2 
using linear finite elements 



Figure 5.6: The cross-section of 
an adaptive mesh of Example 2 
using quadratic finite elements 




Figure 5.7: The convergence 
curves of Example 2 using lin- 
ear finite elements 



Figure 5.8: The convergence 
curves of Example 2 using 
quadratic finite elements 
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where = (0,1) 3 . The exact solution is given by u — sin(7ra:i) sin(7TX2) sin (772:3 

Fig. I5.5l and Fig. I5.6l are two adaptively refined meshes, which show that the 
error indicator is good. It is shown from Fig. 15.71 and Fig. I5.8l that \u — is 
proportional to the a posteriori error estimators, which implies the a posteriori 
error estimators given in section 14.21 are efficient. Besides, similar conclusions 
to that of Example 1 can be obtained from Fig. 15.71 and Fig. 15.81 too. 

Example 3. Consider the Kohn-Sham equation for helium atoms: 



1 



-A 



2 

H 



p(v) 
\x-y\' 



dy + V xc I u = Xu in 



and L 3 I it 1 2 = 1, here p = 2|u| 2 . In our computation of the ground state energy, 
we solve the following nonlinear eigenvalue problem: Find (X,u) e R x Hq(Q) 



such that J n \u\ 2 dx = 1 and 



-A 



2 

H 



p(y) 
\x - y\ 



dy + V xc )u 



Xu in J7, 



u = on dil, 



(5.52) 



where CI = (-10.0, 10. 0) 3 , and V xc {p) = -fa(fp)* with a = 0.77298. Since 
(|5.52j) is a nonlinear eigenvalue problem, we need to linearize and solve them 
iteratively, which is called the self-consistent approach [H [Til [THl [5T] . In our 
computation, a Broyden-type quasi-Newton method |22j were used. 

In 1989, White [27] computed helium atoms over uniform cubic grids and 
obtained ground state energy -2.8522 a.u. by using 500,000 finite clement bases. 
While the ground state energy of helium atoms in Software package fhi98PP [11] 
is -2.8346 a.u., which we take as a reference. 
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Figure 5.9: The ground state en- 
ergy using linear finite elements 



Figure 5.10: The ground state 
energy using quadratic finite ele- 
ments 
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Z=0.0 



Z=0.0 



Figure 5.11: The cross-section of 
an adaptive mesh of Example 3 
using linear finite elements 



Figure 5.12: The cross-section of 
an adaptive mesh of Example 3 
using quadratic finite elements 



Ihe a posteriori estimator 
- a line with slope -1/3 




the a posteriori estimator 
- a line with slope -2/3 




number ol degrees of freedom 



number df degrees ot freeddm 



Figure 5.13: The convergence 
curve of Example 3 using linear 
finite elements 



Figure 5.14: The convergence 
curve of Example 3 using 
quadratic finite elements 
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Our results are displayed in Fig. 15.01 Fig. 15.101 Fig. 15.111 Fig. 15.121 Fig. 
15.131 and Fig. 15.141 It is seen from Fig. 15.101 that the ground state energy in 
our computation is close to the reference with less 100,000 degrees of freedom 
when the quadratic finite element discretization is used. Some cross-sections of 
the adaptively refined meshes are displayed in Fig. I5.11l and Fig |5.12l Since we 
do not have the exact solution, we list the convergence curves of the a posteriori 
error estimators in Fig. 15.131 and Fig. 15.141 only. It is shown from these figures 
that the a posteriori error estimators given in section [4.31 are efficient. 
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